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Abstract. The Helmholtz free energy density is parametrized as a function of temperature 
and baryon density near the chiral critical point of QCD. The parametrization incorporates 
the expected critical exponents and amplitudes. An expansion away from equilibrium states is 
achieved with Landau theory. This is used to calculate the probability that the system is found 
at a density other than the equilibrium one. Such fluctuations are predicted to be very large in 
heavy ion collisions. 



1. Introduction 



> 

o 

00 . If the up and down quark masses are zero and the strange quark mass is not, the transition 

^D I between quark-gluon plasma and hadronic matter may be first or second order at zero baryon 

l/^ ■ chemical potential. If the up and down quark masses are small enough, there may be a phase 

^P i transition at sufficiently large baryon chemical potential. This phase transition is predicted 

to be in the same universality class as liquid-gas phase transitions and the 3D Ising model. 
Diverse studies suggest that there is a curve of first-order phase transition in the /i-T plane that 
terminates in a second-order phase transition at some critical point (^c, Tc). The location of this 
^ ' chiral critical point has been estimated using various effective field theory models, such as the 

Namu Jona-Lasinio model [I]-[3], a composite operator model [1], a random matrix model [S], a 
linear a model [3], an effective potential model [6j, and a hadronic bootstrap model [7], as well 
as various implementations of lattice QCD [8]-[llj. Reviews have been written by Stephanov 
fn\ and Mohanty [T3] . 

This subject is of experimental interest because collisions between heavy nuclei at medium 
to high energy, such as at the future Facility for Antiproton and Ion Research (FAIR), or in 
future low energy runs at the Relativistic Heavy Ion Collider (RHIC), may provide experimental 
information on the phase diagram in the vicinity of a critical point. One characteristic signature 
would be large fiuctuations on an event by event basis |14j-[17j. 

The goal of this work [18] is to understand the equation of state of QCD near the chiral 
critical point and some of its implications for high energy heavy ion collisions. Among the 
requirements are the incorporation of critical exponents and amplitudes and to have sensible 
limits as temperature T — t- and baryon chemical potential /i — > 0. This is accomplished by 
parametrizing the Helmholtz free energy density as a function of T and baryon density n in an 
appropriate fashion. Perhaps the closest work that addressed some of these issues blended a 
parameterization of the 3D Ising model equation of state into quark and hadron equations of 
state |19]. These two parameterizations can perhaps be viewed as alternatives which provide 
some idea as to the range of uncertainty in how to describe matter near the chiral critical point. 



2. Critical curve 

Numerous studies suggest that the curve in the T vs. fi plane separating the quark-gluon and 
hadron phases is approximately quadratic. 
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See figure 1. The parameters Tq = 180 MeV and /^o = 1230 MeV were chosen to match 
theoretical estimates and experimental data for the transition between the two phases at ^ = 
and T = 0, respectively. In this example the critical point is chosen to have the temperature 
Tc = 140 MeV. Along the solid portion of the curve the transition is first order. It terminates at 
a second order transition at the critical point. The dashed portion of the curve indicates that 
there is no true thermodynamic phase transition separating the phases, only a rapid crossover. 
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Figure 1. Temperature versus 
baryon chemical potential from the 
parameterization given in the text. 
The critical point lies somewhere 
along this curve. 



Figure 2. The pressure, entropy 
density, and energy density, nor- 
malized so that they all have the 
same asymptotic value, at /U = 0. 
The parameterization is given in 
the text. 



We will need the values of the pressure Pc, energy density ec, baryon density Uc., and entropy 
density Sc at the critical point. The parametrization of the equation of state near the critical 
point, to be given below, is normalized to these values. Let us assume that the high energy 
density equation of state can be parameterized as 



P 



AiT^ + A2ii'T' + Ao/ - B2T^T' - B^T^ (2) 

The dimensionless coefficients An are adjusted to match a free gas of gluons and 2.5 flavors of 
massless quarks (the 1/2 to approximately take into account the mass of the strange quark). 
The dimensionless coefficients Bn are adjusted to reproduce lattice results near the crossover 
when // = [20]- [22], and to make the pressure a constant along the critical curve when ([1]) is 
inserted. The resulting equation of state for n = and T > Tq is shown in figure 2. 



3. Equation of state near the critical point 

Given the Helmhotz free energy density f{n,T) the other thermodynamic quantites follow 
uniquely from it, as summarized below. 

2d_ n\n,Ty 

dn \ n 



P = n^T^ 



(3) 



. ^ -'-^ (5) 

e = f{n,T) + Ts{n,T) (6) 

These satisfy the thermodynamic identity e = —P + Ts + fin. It is useful to define the 
dimensionless variables t = (T — Tc)/Tc and r] = [n — Uf^jnc which are a measure of distance 
from the critical point in temperature and density. Near the critical point the heat capacity cy 
has the divergent behavior 

_ ds{n,T) f C-{-tY^ whent^O- 
^^ dT 1 c+t-° whent^O+, ^^ 

while the thermal conductivity kt and baryon susceptibility xb have the divergent behavior 
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K_(— t) "^ when t — )■ . . 

K+t-'y when t ^ 0+ . ^^ 



Along the coexistence curve the density difference between the higher density liquid phase 
(quarks and gluons) and the lower density gas phase (hadrons) goes to zero as 

ni-ng^ i-tf. (9) 

Along the critical isotherm the pressure behaves as 

P-Pc^\rifsign{r]). (10) 

Here a, /3, 7 and 5 are critical exponents. These exponents are related by a + 2/3 + 7 = 2 and 
7 = f3{5 — 1). Mean field theories normally give a = 0, /3 = 1/2, 7 = 1, and 6 = 3. Typical 
fluids are measured to have a < 1, /3 Ri 1/3, 1.2 < 7 < 1.3, and 4 < J < 5 [23j. The 3D 
Ising model has a = 0.11, /3 = 0.325, 7 = 1.24, and 5 = 4.815 [24], which are the values used 
here. Furthermore, the critical amplitudes are related by k+/k^ ~ 5 and c+/c- ~ 0.5, which is 
universal to all theories within the same class. See [25]. 

The challenge now is to parametrize /(??,t) so that it yields the correct critical behavior 
near r] = t = and still gives sensible results away from the critical point. Mean field theories 
normally result in a Taylor series expansion in integral powers of rj with coefficients that depend 
on t. The simplest way to incorporate the correct critical behavior is with a modified expansion 
of the form 

/ = hit) + /i(t)r? + f2{tW + Ut)\vr. (11) 

The power a = 6 + 1 = 5.815. Higher powers of 77 could be included but they would not affect 
the critical behavior and so are discarded here for simplicity. The coefficient functions must 
have the following behavior. 

/o(t) = \ ''Dli _ " Va if t > . ^^"^^ 

^(1 + i)'- (13) 

3 

Here fo{t) and /2(i) are smooth functions of t with /2(0) = 0. They are parametrized so that 
/(n, T) has sensible limits at T = and at n = 0. The coefficient f^ is assumed to be constant 
and proportional to P^, for definiteness we shall take fa- = 5Pc ~ 512 MeV/fm^. Results for the 
critical curves, latent heat, thermal conductivity, and heat capacity are shown in figures 3 to 6. 
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Figure 3. The solid curve 

denotes coexistence between high 
and low density phases and the 
dashed curve denotes the limits of 
metastability. The scaled curves 
are independent of the choice of 
critical temperature and density. 
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Figure 4. The latent heat per 
unit volume versus temperature for 
three choices of critical tempera- 
ture. 
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Figure 5. The isothermal com- 
pressibility. For t < they are eval- 
uated along the coexistence curve 
while for t > it is evaluated at 
the critical density. The curves for 
different critical temperature fall on 
top of one another. 
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Figure 6. The heat capacity per 
unit volume. For i < they 
are evaluated along the coexistence 
curve while for i > it is evaluated 
at the critical density. When 
divided by the entropy density at 
the critical point the results are 
nearly independent of the choice of 
critical temperature. 
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4. Fluctuations 

In any system of finite size there will be thermal fluctuations. These fluctuations become large 
as the critical point is approached. Fluctuations can be studied with Landau theory [261 127] . 



The question is: What is the probabihty to find a system of volume V with a baryon density 
different from the equihbrium one? The expansion away from equilibrium states is determined 
by the thermodynamic potential 



nifi,T;r])-no{fi,T)= (f^-ncf^)r] + f2V^ + fM'' V. 



(15) 



Here i^o{fi,T) = (/o — ncfi)V. Along the coexistence curve /i = ndJ'- The relative probability 
to be at a density other than the equilibrium one, along the coexistence curve, is 



^(^)/^(^0 = exp (-Af]/T) 



(16) 



where 



An 



f2[v-vf)+fA\vr 



\m\ 



V. 



(17) 



For purposes of illustration the volume is taken to be 400 fm^. The value of 400 fm^ is about as 
large as one can imagine for high energy nuclear collisions. Considering that the critical density 
is estimated to be about 5no ~ 0.75 baryons/fm^, this would mean that about 300 baryons 
participate in the fluctuation. That is a substantial fraction of the total in Au+AU, Pb+Pb, or 
in U+U collisions. Smaller volumes would have even larger fluctuations. 
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Figure 7. Thermodynamic po- 
tential as a function of r/ near the 
critical point when the volume is 
400 fm^. The stable phases are lo- 
cated at the minima of the poten- 
tial. Four different temperatures 
are shown, with the solid curve rep- 
resenting the critical temperature. 
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Figure 8. The probability to find 
the system at a particular density 
relative to the equilibrium densities 
at phase coexistence. 



The thermodynamic potential is shown in figure 7 for several different temperatures. Phase 
coexistence occurs below T^, the two phases are located at the minima of the potential. They 
are not very deep, indicating that fluctuations do not cost much free energy. The potential is so 
flat because f2{Tc) = and because the power o" « 6 is so large. The shallowness of the potential 
is reflected in the probabilities, as shown in figure 8. The probability to find the system with 
a density anywhere between ±25% of the critical density is greater than 50%. Fluctuations are 
large! They would be even larger for smaller volume systems. 



5. Conclusions 

An equation of state valid in the vicinity of the chiral critical point has been constructed. 
It incorporates correct values of the critical exponents and amplitudes. Since only certain 
properties of the equation of state are universal, there is some freedom to vary the noncritical 
functional dependence on temperature and density. Work on extending the equation of state to 
a wider range of T and // is underway. 

The Landau theory of fluctuations away from equilibrium states was used to determine the 
magnitude of the fluctuations one might expect in heavy ion collisions. These fluctuations are 
quite large, partly due to finite volume effects but mostly because the critical exponent 5 is much 
larger than in mean field theories. This flattens the Landau free energy as a function of density 
away from the equilibrium densities and so decreases the cost to fluctuate away from them. 

It will require careful thought as to how to incorporate fluctuations near the critical point 
into dynamical simulations of heavy ion collisions. What is the appropriate way to describe the 
transition in heavy ion collisions? Is it nucleation [28], spinodal decomposition [29], or something 
else? What are the best experimental observables and can they be measured at RHIC and/or 
FAIR? The future of this topic is exciting! 
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